Vehicular Traffic: A System of Interacting Particles Driven Far Prom Equilibrium* 
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In recent years statistical physicists have developed discrete "particle- hopping" models of vehicular 
traffic, usually formulated in terms of cellular automata, which are similar to the microscopic models 
of interacting charged particles in the presence of an external electric field. Concepts and techniques 
of non-equilibrium statistical mechanics are being used to understand the nature of the steady states 
and fluctuations in these so-called "microscopic" models. In this brief review we explain, primarily 
to the nonexperts, these models and the physical implications of the results. 



I. INTRODUCTION 

Are you surprised to see an article on vehicular traffic 
in this special section of current science where physicists 
are supposed to report on some recent developments in 
the area of dynamics of noncquilibrium statistical sys- 
tems? Aren't civil engineers (or, more specifically, traffic 
engineers) expected to work on traffic? Solving traffic 
problems would become easier if one knows the funda- 
mental laws governing traffic flow and traffic jam. For 
almost half a century physicists have been trying to de- 
velope a theoretical framework of traffic science extend- 
ing concepts and techniques of statistical physics [H-p|. 
The main aim of this brief review is to show how these 
attempts, particularly the recent ones, have led to deep 
insight in this frontier area of inter-disciplinary research. 

The dynamical phases of systems driven far from equi- 
librium are counterparts of the stable phases of systems 
in equilibrium. Let us first pose some of the questions 
that statistical physicists have been addressing in or- 
der to discover the fundamental laws governing vehicular 
traffic. For example, 

(i) What are the various dynamical phases of traffic? 
Does traffic exhibit phase-coexistence, phase transition, 
criticality or self-organized criticality and, if so, under 
which circumstances? 

(ii) What is the nature of fluctuations around the steady- 
states of traffic? 

(iii) If the initial state is far from a stationary state of 
the driven system, how does it evolve with time to reach 
a truely steady-state? 

(iv) What are the effects of quenched disorder (i.e., time- 
independent disorder) on the answers to the questions 
posed in (i)-(iii) above? 

The microscopic models of vehicular traffic can find prac- 
tical applications in on-line traffic control systems as well 
as in the planning and design of transportation network. 

There are two different conceptual frameworks for 
modellig vehicular traffic. In the "coarse-grained" fluid- 
dynamical description, the traffic is viewed as a com- 
pressible fluid formed by the vehicles but these individual 
vehicles do not appear explicitly in the theory. In con- 
trast, in the "microscopic" models traffic is treated as a 



system of interacting "particles" driven far from equilib- 
rium where attention is explicitly focussed on individual 
vehicles each of which is represented by a "particle"; the 
nature of the interactions among these particles is de- 
termined by the way the vehicles influence each others' 
movement. Unlike the particles in a gas, a driver is an 
intelligent agent who can "think", make individual deci- 
sions and "learn" from experience. Nevertheless, many 
general phenomena in traffic can be explained in general 
terms with these models provided the behavioural effects 
of the drivers are captured by only a few phenomenolog- 
ical parameters. 

The conceptual basis of the older theoretical ap- 
proaches are explained briefly in section II. Most of the 
"microscopic" models developed in the recent years are 
" particle- hopping" models which are usually formulated 
using the language of cellular automata (CA) j^. The 
Nagel-Schreckenberg (NaSch) g model and the Biham- 
Middleton-Levine (BML) M model, which are the most 
popular CA models of traffic on idealized highways and 
cities, respectively, have been extended by several au- 
thors to develope more realistic models. Some of the 
most interesting aspects of these recent developments are 
discussed in the long sections III and IV. The similarities 
between various particle-hopping models of traffic and 
some other models of systems, which are also far from 
equilibrium, are pointed out in section V followed by the 
concluding section VI. 



II. OLDER THEORIES OF VEHICULAR 
TRAFFIC 

A. Fluid-dynamical Theories of vehicular traffic 

In traffic engineering, the fundamental diagram depicts 
the relation between density c and the flux J, which 
is defined as the number of vehicles crossing a detector 
site per unit time fM]. Because of the conservation of 
vehicles, the local density c{x;t) and local flux J{x;t) 
satisfy the equation of continuity which is the analogue 
of the equation of continuity in the hydrodynamic the- 
ories of fluids. In the early works ||l2] it was assumed 



(i) that the flux (or, equivalently, the velocity) is a func- 
tion of the density and (ii) that, following any change 
in the local density, the local speed instantaneously re- 
laxes to a magnitude consistent with the new density at 
the same location. However, for a more realistic descrip- 
tion of traffic, in the recent fluid-dynamical treatments 
||l3|-p^ of traffic an additional equation (the analogue 
of the Navier-Stokes equation for fluids), which describes 
the time-dependence of the velocity V{x; t), has been con- 
sidered. This approach, however, has its limitations; for 
example, viscosity of trafhc is not a directly measureable 
quantity. 



B. Kinetic theory of vehicular traffic 

In the kinetic theory of traffic, one begins with the 
basic quantity g(x,v,w;t) dx dv dw which is the num- 
ber of vehicles, at time i, located between x and x + dx, 
having actual velocity between v and v + dv and desired 
velocity between w and w + dw. In this approach, the 
fundamental dynamical equation is the analogue of the 
Boltzmann equation in the kinetic theory of gases 0. 
Assuming reasonable forms of "relaxation" and "interac- 
tion" , the problem of traffic is reduced to that of solving 
the Boltzmann-like equation, a formidable task, indeed 



C. Car- following theories of vehicular traffic 

In the car-following theories one writes, for each indi- 
vidual vehicle, an equation of motion which is the ana- 
logue of the Newton's equation for each individual parti- 
cle in a system of interacting classical particles. In New- 
tonian mechanics, the acceleration may be regarded as 
the response of the particle to the stimulus it receives in 
the form of force which includes both the external force as 
well as those arising from its interaction with all the other 
particles in the system. Therefore, the basic philisophy 
of the car-following theories ||l|,y can be summarized by 
the equation 



[Response\n oc [Stimulus]n 



(1) 



for the n-th vehicle (n = 1,2, ...). The constant of pro- 
portionality in the equation (ni) can be interpreted as a 
measure of the sensitivity coefficient of the driver; it indi- 
cates how strongly the driver responds to unit stimulus. 
Each driver can respond to the surrounding traffic con- 
ditions only by accelerating or decelerating the vehicle. 
The stimulus and the sensitivity factor are assumed to be 
functions of the position and speed of the vehicle under 
consideration and those of its leading vehicle. Different 
forms of the equations of motion of the vehicles in the dif- 
ferent versions of the car-following models arise from the 
differences in their postulates regarding the nature of the 



stimulus. In general, the dynamical equations for the ve- 
hicles in the car-following theories are coupled non- linear 
differential equations p9^E5[ and thus, in this "micro- 
scopic" approach, the problem of traffic flow reduces to 
problems of nonlinear dynamics. 



III. CELLULAR-AUTOMATA MODELS OF 
HIGHWAY-TRAFFIC 

In the car-following models space is treated as a con- 
tinuum and time is represented by a continuous variable 
t while velocities and accelerations of the vehicles are 
also real variables. However, most often, for numerical 
manipulations of the differential equations of the car- 
following models, one needs to discretize the continuous 
variables with appropriately chosen grids. In contrast, in 
the CA models of traffic not only time but also the po- 
sition, speed, and acceleration of the vehicles are treated 
as discrete variables. In this approach, a lane is repre- 
sented by a one-dimensional lattice. Each of the lattice 
sites represents a "cell" which can be either empty or 
occupied by at most one "vehicle" at a given instant of 
time (see fig.l). 



FIG. 1. A typical configuration in a CA model. The 
number in the upper right corner is the speed of the vehicle. 

At each discrete time step t —f t + 1, the state of the 
system is updated following a well defined prescription. 



A. The Nagel-Schreckenberg model of highway 
traffic: 

In the NaSch model, the speed V of each vehicle can 
take one of the Vmax + 1 allowed integer values V — 
0,1, ...,Vmax- Suppose, A„ and Vn denote the position 
and speed, respectively, of the n-th vehicle. Then, dn = 
ATn+i — Xn, is the gap in between the n-th vehicle and 
the vehicle in front of it at time t. At each time step 
t ^ t + 1, the arrangement of the N vehicles on a finite 
lattice of length L is updated in parallel according to the 
following " rules" : 

Step 1: Acceleration. If, Vn < Vmax, the speed of the n- 
th vehicle is increased by one, but Vn remains unaltered 

if V„ = Vmax, i.e., Vn ^ TOZn(K -I- 1, Vmax)- 

Step 2: Deceleration (due to other vehicles). If dn < Vn, 
the speed of the n-th vehicle is reduced to (i„ — 1, i.e., 
Vn -+ min{Vn,dn - 1). 

Step 3: Randomization. If Vn > 0, the speed of the n-th 
vehicle is decreased randomly by unity with probability 



p but Vn does not change if Vn — 0, i.e., Vn -^ max{Vn — 
1, 0) with probabihty p. 

Step 4- Vehicle movement. Each vehicle is moved for- 
ward so that Xn — > Xn + Vn- 

The NaSch model is a minimal model in the sense that 
all the four steps are necessary to reproduce the basic fea- 
tures of real traffic; however, additional rules need to be 
formulated to capture more complex situations. The step 
1 reflects the general tendency of the drivers to drive as 
fast as possible, if allowed to do so, without crossing the 
maximum speed limit. The step 2 is intended to avoid 
collision between the vehicles. The randomization in step 
3 takes into account the different behavioural patterns of 
the individual drivers, especially, nondeterministic accel- 
eration as well as overreaction while slowing down; this 
is crucially important for the spontaneous formation of 
traffic jams. So long as p ^ 0, the NaSch model may be 
regarded as stochastic C A . For a realistic description 
of highway traffic |8| , the typical length of each cell should 
be about 7.5m and each time step should correspond to 
approximately 1 sec of real time when Vmax = 5- 

The update scheme of the NaSch model is illustrated 
with a simple example in fig.0. 

Configuration at time t: 



2 




1 






1 








a) Acceleration: 
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b) Braking: 
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c) Randomization (p=l/3): 
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d) Driving (= configuration at time t + 1): 



c^ 



c^ 



FIG. 2. Step-by-step example for the application of the 
update rules. We have assumed Vmax = 2 and p = 1/3. 
Therefore on average one third of the cars qualifying will slow 
down in the randomization step. 




V 



FIG. 3. Typical spacc-tiinc diagrams of the NaSch model 
with Vmax = 5 and (a) p = 0.25, c = 0.20, (b) p = 0.0, c = 0.5. 
Each horizontal row of dots represents the instantaneous posi- 
tions of the vehicles moving towards right while the successive 
rows of dots represent the positions of the same vehicles at 
the successive time steps. 

Space-time diagrams showing the time evolutions of 
the NaSch model demonstrate that no jam is present at 
sufficiently low densities, but spontaneous fluctuations 
give rise to traffic jams at higher densities (fig. 3(a)). 
From the fig. 3(b) it should be obvious that the intrinsic 
stochasticity of the dynamics |^] , arising from non-zero p, 
is essential for triggering the jams |P,|26[|. 

The use of parallel dynamics is also important. In con- 
trast to a random sequential update, it can lead to a 
chain of overreactions. Suppose, a vehicle slows down 
due the randomization step. If the density of vehicles is 
large enough this might force the following vehicle also to 
brake in the deceleration step. In addition, if p is not too 
small, it might brake even further in Step 3. Eventually 
this can lead to the stopping of a vehicle, thus creating 
a jam. This mechanism of spontaneous jam formation is 
rather realistic and cannot be modelled by the random 
sequential update. 



B. Relation betw^een the NaSch model and ASEP 

In the NaSch model with Vmax — 1 every vehicle moves 
forward with probability q = I ~ p in the time step t + 1 
if the site immediately in front of it were empty at the 
time step t; this, is similar to the fully asymmetric sim- 
ple exclusion process (ASEP) p7|-p9[ where a randomly 
chosen particle can move forward with probablity q if 
the site immediately in front is empty. But, updating is 
done in parallel in the NaSch model whereas that in the 
ASEP is done in a random sequential manner. Neverthe- 
less, the special case of Vmax — 1 for the NaSch model 
achieves special importance from the fact that so far it 
has been possible to derive exact analytical results for 
the NaSch model only in the special limits (a) Vnax = 1 
and arbitrary p and (b) p = and arbitrary Vmax ■ 



C. NaSch model in the deterministic limits 

If p = 0, the system can self-organize so that at low 
densities every vehicle can move with Vmax and the cor- 
responding flux is cVmax] this is, however, possible only 
if enough empty cells are available in front of every ve- 
hicle, i.e., for c < c^* = {Vmax + 1)"^ and the cor- 
responding maximum flux is J^^^ = Vmax /{Vmax + !)• 
On the other hand, for c > d^*' , the flow is hmited by 
the density of holes. Hence, the fundamental diagram 
in the deterministic limit p = of the NaSch model 
(for any arbitrary Vmax) is given by the exact expression 
J = min{cVmax, {I- c)). 

Aren't the properties of the NaSch model with max- 
imum allowed speed Vmax, in the deterministic limit 
p = 1, exactly identical to those of the same model with 
maximum allowed speed Vmax — 1? The answer to the 
question posed above is: NO; iip = 1, all random initial 
states lead to J = in the stationay state of the NaSch 
model irrespective of Vmax and c! 



D. Analytical Theory for the NaSch Model 



For example, the configuration shown in fig. 4 is a GoE 
stateQ because it could occur at time t only if the two 
vehicles occupied the same cell simultaneously at time 
t — 1. The naive mean- field theory mentioned above does 
not exclude the GoE states. The exact expression, given 
in the next subsection, for the flux in the steady-state 
of the NaSch model with Vmax — 1 can be derived by 
merely excluding these states from consideration in the 
naive mean- field theory ||3^, thereby indicating that the 
only source of correlation in this case is the parallel up- 
dating. But, for Vmax > 1, there are other sources of 
correlation because of which exclusion of the GoE states 
merely improves the naive mean- field estimate of the flux 
but does not yield exact results | |3l| ]. 

A systematic improvement of the naive mean-field the- 
ory of the NaSch model has been achieved by incorporat- 
ing short-ranged correlations through cluster approxima- 
tions. We define a n-cluster to be a collection of n suc- 
cessive sites. In the general n-cluster approximation, one 
divides the lattice into "clusters" of length n such that 
two neighbouring clusters have n—1 sites in common (see 
fig.5). 



In the "site-oriented" theories one describes the state 
of the finite system of length L by completely specifying 
the state of each site. In contrast, in the "car-oriented" 
theories the state of the traffic system is described by 
specifying the positions and speeds of all the N vehicles 
in the system. In the naive mean-field approximation one 
treats the probabilities of occupation of the lattice sites 
as independent of each other. In this approximation, for 
example, the steady-state flux for the NaSch model with 
Vmax = 1 and periodic boundary conditions, one gets pG] 



J = qc{l — c) 



(2) 



It turns out [ pO[ that the naive mean-field theory un- 
derestimates the flux for all Vmax- Curiously, if instead 
of parallel updating one uses the random sequential up- 
dating, the NaSch model with Vmax = 1 reduces to the 
ASEP for which the equation (0) is known to be the exact 
expression for the corresponding flux (see, e.g., [gj)! 

What are the reasons for these differences arising 
from parallel updating and random sequential updating? 
There are "garden of Eden" (GoE) states (dynamically 
forbidden states) [pl| of the NaSch model which cannot 
be reached by the parallel updating whereas no state is 
dynamically forbidden if the updating is done in a ran- 
dom sequential manner. 



b) 




c) 




FIG. 5. Decomposition of a lattice into (a) 1-clusters (b) 
2-clusters and (c) 3-clusters in the cluster-theoretic approach 
to the NaSch model. 

If n = 1, then the 1-cluster approximation can be re- 
garded as the naive mean-field approximation. You can 
easily verify, for example, in the special case of Vmax — 1, 
that the state of the 2-cluster at time t+1 depends on the 
state of the 4-cluster at time t, which, in turn, depends 
on the state of a larger cluster at time t — 1 and, so on. 



FIG. 4. A GoE state for the NaSch model with V„iax > 2. 



The configuration shown in fig.l is also a GoE state! 



Therefore, one needs to make an approximation to trun- 
cate this hierarchy in a sensible manner. For example, 
in the 2-cluster approximation for the NaSch model with 
Vmax = 1, the 4-cluster probabilities are approximated 
in terms of an appropriate product of 2-cluster probabil- 
ities. Thus, in the n-cluster approximation ||30|] a cluster 
of n neighbouring cells are treated exactly and the cluster 
is coupled to the rest of the system in a self-consistent 
way. 

Carrying out the 2-cluster calculation |3^ for Vmax — 1 
one not only finds an effective particle-hole attraction 
(particle-particle repulsion), but also obtains the exact 
result 



Jic,p)^^[l-^l-Aqc{l-. 



(3) 



for the corresponding flux. But one gets only ap- 
proximate results from the 2-cluster calculations for all 
Vmax > 1 (see ||32J for higher order cluster calculations 
for Vmax — 2 and comparison with computer simulation 
data) . 
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FIG. 6. The fundamental diagram in the Nasch model for 
(a) V,nax = 1 and (b) V^nax > 1, both for p = 0.25. The 
data for all Vmax > 1 have been obtained through computer 
simulations. 

Let us explain the physical origin of the generic shape 
of the fundamental diagrams shown in fig. 6. At suf- 
ficiently low density of vehicles, practically "free flow" 
takes place whereas at higher densities traffic becomes 
"congested" and traffic jams occur. So long as c is suffi- 
ciently small, the average speed {V) is practically inde- 
pendent of c as the vehicles are too far apart to interact 
mutually. However, a faster monotonic decrease of (V) 
with increasing c takes place when the forward movement 
of the vehicles is strongly hindred by others because of 
the reduction in the average separation between them. 
Because of this trend of variation of (V) with c, the flux 
J = (cV) exhibits a maximum pXJ ] at c.,„; for c < Cm, in- 
creasing c leads to increasing J whereas for c > Cm sharp 
decrease of (V) with increase of c leads to the overall 
decrease of J. 

An interesting feature of the expression (Q) is that the 
flux is invariant under charge conjugation, i.e., under the 
operation c -^ (1 — c) which interchanges particles and 
holes. Therefore, the fundamental diagram is symmetric 
about c = 1/2 when Vmax = 1 (see fig. 6a). Although 
this symmetry breaks down for all Vmax > 1 (see fig. 6b), 
the corresponding fundamental diagrams appear more re- 
alistic. Moreover, for given p, the magnitude of c,„ de- 



creases with increasing Vmax as the higher is the Vmax the 
longer is the effective range of interaction of the vehicles 
(see fig. 6b). Furthermore, for Vmax = 1, dux merely de- 
creases with increasing p (see eqn(|3|)), but remains sym- 
metric about c — 1/2 — Cm- On the other hand, for all 
Vmax > 1 , increasing p not only leads to smaller fiux but 
also lowers Cm- 



E. Spatio-temporal organization of vehicles 

The distance from a selected point on the lead vehicle 
to the same point on the following vehicle is defined as 
the distance- headway (DH) p(| . In order to get infor- 
mation on the spatial organization of the vehicles, one 
can calculate the DH distribution Vdhi^X) by following 
either a site-oriented approach |34[| or a car-oriented ap- 
proach [p3| if AXj = Xj — Xj-i, i.e., if the number of 
empty lattice sites in front of the j-th vehicle is identified 
as the corresponding DH. At moderately high densities, 
Vdhi'^^X) exhibits two peaks; the peak at AX = 1 is 
caused by the jammed vehicles while that at a larger AX 
corresponds to the most probable DH in the free-flowing 
regions. 

The time-headway is defined as the time interval be- 
tween the departures (or arrivals) of two successive vehi- 
cles recorded by a detector placed at a fixed position on 
the highway |10(| . The time- headway distribution con- 
tains information on the temporal organization. Sup- 
pose, Vm{ti) is the probability that the following vehicle 
takes time ti to reach the detector, moving from its ini- 
tial position where it was located when the leading vehicle 
just left the detector site. Suppose, after reaching the de- 
tector site, the following vehicle waits there for r—ti time 
steps, either because of the presence of another vehicle 
in front of it or because of its own random braking; the 
probability for this event is denoted by Q{t ~ ti \ti). The 
distribution Vth{T), of the time-headway t, can be ob- 
tained from [||§] Vthir) = Eli=\ ^m(*i)Q(^ " ^il^i)- 
The most-probable time-headway, when plotted against 
the density, exhibits a minimum p6[ ; this is consistent 
with the well known exact relation J — l/Tav between 
fiux and the average time- headway, Tav 

Is there a phase transition from " free- flowing" to " con- 
gested" dynamical phase of the NaSch model? No satisfa- 
cory order parameter has been found so far fS^JS^ , except 
in the deterministic limit p9] . The possibility of the exis- 
tence of any critical density in the NaSch model is ruled 
out by the observations |^,^,^,Q that, for all non- 
zero p, (a) the equal-time correlation function decays ex- 
ponentially with separation, and (b) the relaxation time 
and lifetimes of the jams remain finite. This minimal 
model of highway traffic also does not exhibit any first 
order phase transition and two-phase co-existence Bm. 



F. Extensions of the NaSch model and practical 
applications 



IV. CELLULAR-AUTOMATA MODELS OF 
CITY-TRAFFIC 



In recent years some other minimal models of traffic on 
highways have been developed by modifying the updating 
rules of the NaSch model fpl-Bsj . In the cruise control 
limit of the NaSch model |44[ | the randomization step is 
applied only to vehicles which have a velocity V < Vmax 
after step 2 of the update rule. Vehicles moving with 
their desired velocity Vmax are not subject to fluctua- 
tions. This is exactly the effect of a cruise-control which 
automatically keeps the velocity constant at a desired 
value. Interestingly, the cruise-control limit of the NaSch 
model exhibits self-organized criticality ^.Q. Besides, 
a continuum limit of the NaSch model has also been con- 
sidered (I?). 

The vehicles which come to a stop because of hindrance 
from the leading vehicle may not be able to start as soon 
as the leading vehicle moves out of its way; it may start 
with a probability qs <!■ When such possibilities are in- 
corporated in the NaSch model, the " slow-to-start" rules 
[p8|-p2[ can give rise to metastable states of very high 
flux and hysteresis effects as well as phase separation of 
the traffic into a "free-flowing" phase and a "mega-jam". 

The bottleneck created by quenched disorder of the 
highway usually slows down traffic and can give rise to 
jams ||5^ , P5|| and phase segregation ||5^ , |55|| . However, a 
different type of quenched disorder, introduced by assign- 
ing randomly different braking probabilities p to different 
drivers in the NaSch model, can have more dramatic ef- 
fects [ p6|j57| ] which are reminiscent of " Bose-Einstein-like 
condensation" in the fully ASEP where particle-hopping 
rates are quenched random variables |pq,p9[. In such 
" Bose-Einstein-like condensed" state finite fraction of the 
empty sites are " condensed" in front of the slowest vehi- 
cle (i.e., the driver with highest p). 

Several attempts have been made to generalize the 
NaSch model to describe traffic on multi-lane highways 
and to simulate traffic on real networks in and around 
several cities pG] . For planning and design of the trans- 
portation network pq] , for example, in a metropolitan 
area [pl|-|63[, one needs much more than just micro- 
simulation of how vehicles move on a linear or square 
lattice under a specified set of vehicle-vehicle and road- 
vehicle interactions. For such a simulation, to begin with, 
one needs to specify the roads (including the number of 
lanes, ramps, bottlenecks, etc.) and their intersections. 
Then, times and places of the activities, e.g., working, 
shopping, etc., of individual drivers are planned. Micro- 
simulations are carried out for all possible different routes 
to execute these plans; the results give informations on 
the efficiency of the different routes and these informa- 
tions are utilized in the designing of the transportation 
network p3-pq]. Some socio-economic questions as well 
as questions on the environmental impacts of the planned 
transportation infrastructure also need to be addressed 
during such planning and design. 



A. The Biham-Middleton-Levin model of city traffic 
and its generalizations 

In the BML model 0, each of the sites of a square 
lattice represent the crossing of a east-west street and 
a north-south street. All the streets parallel to the X- 
direction of a Cartesian coordinate system are assumed 
to allow only single-lane east-bound traffic while all those 
parallel to the y-direction allow only single-lane north- 
bound traffic. In the initial state of the system, vehicles 
are randomly distributed among the streets. The states 
of east-bound vehicles are updated in parallel at every 
odd discrete time step whereas those of the north-bound 
vehicles are updated in parallel at every even discrete 
time step following a rule which is a simple extension of 
the fully ASEP: a vehicle moves forward by one lattice 
spacing if and only if the site in front is empty, otherwise 
the vehicle does not move at that time step. 

Computer simulations demonstrate that a first order 
phase transition takes place in the BML model at a finite 
non- vanishing density c, , where the average velocity of 
the vehicles vanishes discontinuously signalling complete 
jamming; this jamming arises from the mutual blocking 
of the flows of east-bound and north-bound traffic at var- 
ious different crossings |^,Q. Note that the dynamics 
of the BML model is fully deterministic and the random- 
ness arises only from the random initial conditions [ p9| . 

As usual, in the naive mean-field approximation one 
neglects the correlations between the occupations of dif- 
ferent sites 1 70 . However, if you are not interested in 



detailed information on the "structure" of the dynami- 
cal phases, you can get a mean-field estimate of c* by 
carrying out a back-of-the-envelope calculation [[7l|-f73[. 



Vy = V, 



In the symmetric case Cx — Cy, for which Vx ^ Vy = 
c = c^~ 0.343. 

The BML model has been extended to take into ac- 
count the effects of (i) asymmetric distribution of the 
vehicles [ [71|, i.e., Cx 7^ Cy, (ii) overpasses or two- level 
crossings [[72[ that are represented by specifically iden- 
tified sites each of which can accomodate upto a max- 
imum of two vehicles simultaneously, (iii) faulty traffic 
lights [Q (iv) static hindrances or road blocks or vehicles 
crashed in traffic accident, i.e., stagnant points | |75| , [76[ , 
(v) stagnant street where the local density Cg of the ve- 
hicles is initially higher than that in the other streets 
| |77[ (vi) jam- avoiding drive |^^ of vehicles to a neigh- 
bouring street, parallel to the original direction, to avoid 
getting blocked by other vehicles in front, (vii) turning 
of the vehicles from east-bound (north-bound) to north- 
bound (east-bound) streets [^. (viu) a single north- 
bound street cutting across east-bound streets [Q (ix) 
more realistic description of junctions of perpendicular 
streets |8ll|8|, (x) green- waves H. 



B. Marriage of NaSch and BML models 

At first sight the BML model may appear very unreal- 
istic because the vehicles seem to hop from one crossing 
to the next. However, it may not appear so unrealistic if 
each unit of discrete time interval in the BML model is 
interpreted as the time for which the traffic lights remain 
green (or red) before switching red (or green) simultane- 
ously in a synchronized manner, and over that time scale 
each vehicle, which faces a green signal, gets an opportu- 
nity to move from j'-th crossing to the j + 1-th (or, more 
generally Q], to the j + r-th where r > 1). 

However, if one wants to develope a more detailed 
"fine-grained" description then one must first decorate 
each bond [g5[ with Z? — 1 {D > 1) sites to represent D — 1 
cells in between each pair of successive crossings thereby 
modelling each segment of the streets in between succes- 
sive crossings in the same manner in which the entire 
highway is modelled in the NaSch model. Then, one can 
follow the prescriptions of the NaSch model for describ- 
ing the positions, speeds and accelerations of the vehicles 
^ , p6| as well as for taking into account the interactions 
among the vehicles moving along the same street. More- 
over, one should flip the color of the signal periodically 
at regular interval of T (T >> 1) time steps where, dur- 
ing each unit of the discrete time interval every vehicle 
facing green signal should get an opportunity to move 
forward from one cell to the next. Such a CA model of 
traffic in cities has, indeed, been proposed very recently 
| p7[ where the rules of updating have been formulated in 
such a way that, (a) a vehicle approaching a crossing can 
keep moving, even when the signal is red, until it reaches 
a site immediately in front of which there is cither a halt- 
ing vehicle or a crossing; and (b) no grid-locking would 
occur in the absence of random braking. 

A phase transition from the "free-flowing" dynamical 
phase to the completely "jammed" phase has been ob- 
served in this model at a vehicle density which depends 
on D and T. The intrinsic stochasticity of the dynamics, 
which triggers the onset of jamming, is similar to that 
in the NaSch model, while the phenomenon of complete 
jamming through self-organization as well as the final 
jammed configurations (fig. 7) are similar to those in the 
BML model. This model also provides a reasonable time- 
dependence of the average speeds of the vehicles in the 
"free-fiowing" phase [p7|. 
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FIG. 7. A typical jammed configuration of tlie vehicles. 
The east-bound and north-bound vehicles are represented by 
the symbols — » and j, respectively. (A'^ = 5,D = 8) 



synthesis has been pointed out [p8[. Another driven- 
dissipative system, which is also receiving wide atten- 
tion of physicists in recent years, is the granular material 
flowing through a pipe ^,g|. There are some superficial 
similarities between the clustering of vehicles on a high- 
way and particle-particle (and particle-cluster) aggrega- 
tion process [g^l ■ 

The NaSch model with Vmax = 1 can be mapped onto 
stochastic growth models of one-dimensional surfaces in 
a two-dimensional medium. Particle (hole) movement 
to the right (left) correspond to local forward growth 
of the surface via particle deposition. In this scenario a 
particle evaporation would correspond to a particle (hole) 
movement to the left (right) which is not allowed in the 
NaSch model. It is worth pointing out that any quenched 
disorder in the rate of hopping between two adjacent sites 
would correspond to columnar quenched disorder in the 
growth rate for the surface |^ . 

Inspired by the recent success in theoretical studies 
of traffic, some studies of information traffic on the 
computer network (internet) have also been carried out 

m. 



V. RELATION WITH OTHER SYSTEMS AND 
PHENOMENA 

You must have noticed in the earlier sections that some 
of the models of traffic are non-trivial generalizations 
or extensions of the ASEP, the simplest of the driven- 
dissipative systems which are of current interest in non- 
equilibrium statistical mechanics [pSj . Some similarities 
between these systems and a dynamical model of protein 



VI. SUMMARY AND CONCLUSION: 

Nowadays the tools of statistical mechanics are increas- 
ingly being used to study self-organization and emergent 
collective behaviour of complex systems many of which, 
including vehicular traffic, fall outside the traditional do- 
main of physical systems. However, as we have shown in 
this article, a strong theoretical foundation of traffic sci- 
ence can be built on the basic principles of statistical me- 
chanics. In this brief review we have focussed attention 



mainly on the progress made in the recent years using 
" particle- hopping" models, formulated in terms of cel- 
lular automata, and compared these with several other 
similar systems. 
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